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Abstract 

On 

Fhe dual Dyson series [M.Frasca, Phys. Rev. A 58, 3439 (1998)], is used to 
develop a general perturbative method for the study of atom-field interaction 
in quantum optics. In fact, both Dyson series and its dual, through renormal- 
ization group methods to remove secular terms from the perturbation series, 
give the opportunity of a full study of the solution of the Schrodinger equa- 
ls) ■ tion in different ranges of the parameters of the given hamiltonian. In view 
of recent experiments with strong laser fields, this approach seems well-suited 
, to give a clarification and an improvement of the applications of the dressed 
states as currently done through the eigenstates of the atom-field interaction, 
showing that these are just the leading order of the dual Dyson series when 
OO ■ the Hamiltonian is expressed in the interaction picture. In order to exploit 
the method at the best, a study is accomplished of the well-known Jaynes- 



Cummings model in the rotating wave approximation, whose exact solution 
is known, comparing the perturbative solutions obtained by the Dyson series 



£h ' and its dual with the same approximations obtained by Taylor expanding 

the exact solution. Finally, a full perturbative study of high-order harmonic 
generation is given obtaining, through analytical expressions, a clear account 
of the power spectrum using a two-level model, even if the method can be 
successfully applied to a more general model that can account for ionization 
5h ' too. The analysis shows that to account for the power spectrum it is needed 

to go to first order in the perturbative analysis. The spectrum obtained gives 
a way to measure experimentally the shift of the energy levels of the atom 
interacting with the laser field by looking at the shifting of hyper-Raman 
lines. 

PACS: 42.50.Ct, 42.50.Hz, 42.65.Ky, 32.80.-t 
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I. INTRODUCTION 



Recent experiments on atoms using strong laser fields M have shown the 
appearance of a wealth of new effects, e.g. high-order harmonics generation, 
in the interaction between light and atoms. This situation forced researchers 
to find different approaches to describe the outcomes of those experiments. 
Numerical studies of the time-dependent Schrodinger equation || have shown 
that the two-level model still proves to be very useful to describe all the fea- 
tures of harmonics generation ||, even if the rotating wave approximation 
must be abandoned. Indeed, recent work [Qj5| indicates, by comparing re- 
sults from a two-level model using Floquet states and numerical work on the 
Schrodinger equation, that the simple two-level model is fairly effective in de- 
scribing the physical situation at hand. So far, no perturbative solution seems 
to be known of this two-level model beyond Floquet states for the case of a 
strong laser field. But, a study by Meystre of an atom in a Fabry- Perot cavity 
H used the same model of Ref. || and gave a first perturbative analytical 
solution to such a model in a strong coupling regime. In fact, the analytical 
solution given by Meystre and its higher order corrections has been success- 
fully obtained in Ref. |7j, showing that the levels of the atom undergoes a 
shift. Being the same model, now we have at hand a way to observe exper- 
imentally such a shift through hyper-Raman lines in harmonic generation, if 
one is able to properly account for the spectrum. 

An understanding of interaction between an atom and a strong electro- 
magnetic field has been possible in recent years through the introduction of 
the dressed-atom picture [|| . This approach assumes that the field couples the 
levels of the atom in such a way that the interaction is between this "dressed" 
atom and the field itself. The computation of the corresponding dressed 
states, as currently found in literature, involves or the computation of the 
eigenstates and the eigenvalues of the term of interaction between the atom 
and the field in the Hamiltonian either the computation of the eigenstates of 
the full Hamiltonian, taking in this way into account the field too. From a 
physical standpoint the dressed-atom picture is quite general as it assumes 
that the photons of the field surround the atom as to modify the way the atom 
itself responds to the field, then it should concern a fully second quantized 
theory. But, the computation of the eigenstates of the full Hamiltonian or just 
the atom-field interaction term, that we take to be the dressed states, often 
reveals itself as an approximation scheme whose understanding is the main 
aim of this paper. So far, no reason has been known for the nice working of 
such dressed states in applied mathematics. A recently devised approach ||, 
the dual Dyson perturbation series, turns out to be both an explanation and 
an improvement of the computation of dressed states permitting the compu- 
tation of higher order corrections to a leading order solution obtained through 
such dressed states. As a by-product one has a clear physical understanding 
of what are the parameters involved in such approximate dressed states and 
what is going to neglect. So, by this improvement of the computation of 
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dressed states, we are able to find an analytical perturbative solution to the 
two-level model to analyse high-order harmonic generation showing that this 
is a first order effect, that is, the leading order solution found by Meystre 
is not enough to get the right spectrum. Then, the result properly accounts 
for the relevance of population distribution as discussed in Ref. |5] and an 
analytical closed expression is given. 

The dual Dyson series that accounts for the dressed states as defined above 
can be derived from the time-dependent Schrodinger equation by using the 
duality principle in perturbation theory and the quantum adiabatic approxi- 
mation H . In this way one realizes that the dual Dyson series is the same one 
of Ref. [ 10 1 . The results one gets from what should work just for quantum 
adiabatic processes can appear somewhat unexpected, as it will be shown for 
the Jaynes-Cummings model in the rotating wave approximation (RWA) for 
whom an exact solution is known. But, this just agrees with the results of 
Ref. §. 

So, the existence of a dual Dyson series can improve the study of atom- 
field interaction. In fact, one can accomplish a perturbative analysis of mod- 
els in quantum optics in different regions of the parameter space that for a 
Jaynes-Cummings model can be easily identified, when spontaneous emission 
is neglected, with the ratio between the detuning and the Rabi frequency. 
Then, by generalizing the computation of dressed states through the dual 
Dyson series on one side and by the standard Dyson series on the other, we 
can reach the main aim of this paper: A general perturbative method to study 
atom-field interaction in quantum optics at different values of the parameters 
of the Hamiltonian. 

The completeness of our approach is strongly tied with the recent results 
obtained in quantum optics through the renormalization group methods for 
perturbation theory Q. These methods permit the resummation of the so 
called secularities that appear in perturbation theory. Indeed, we are able 
to derive a energy level shift of the atom in high-order harmonic generation 
that has effect on hyper-Raman lines. As shown in Ref. ||, when the two 
levels of the atom are equally populated, only hyper-Raman lines should be 
observed. Then, in view of this situation, such an energy level shift turns out 
to be significant. 

It should be pointed out that, although the extension of this approach 
to the method of the master equation M should be straigthforward, it is not 
considered in this paper. So, e.g. the effect of vacuum fluctuations of the field 
modes is neglected. 

The paper is so structured. In sec. II we give a general description of the 
methods and show why the eigenstates of the perturbation are important for 
strong fields. In sec. Ill a study of the Jaynes-Cummings model in RWA is 
accomplished in order to have a pedagogical description of the methods and 
a comparation with an exact solution. In sec. IV the question of high-order 
harmonic generation is discussed through the methods so far introduced. 
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II. A GENERAL METHOD FOR PERTURB ATIVE ANALYSIS 



A. General theory 

In Ref. H we have introduced the duality principle in perturbation theory. 
By duality we mean that, for a given differential equation, it is possible to 
compute both a perturbation series in A and j, being A the characteristic 
parameter of the equation. This is accomplished by a proper choice of the 
leading order equation. So, e.g. for the Duffing equation 

x + x + Xx 3 = 

one can compute a series in A and j by taking, at leading order, in the former 
case x + x = and in the latter case x + Ax 3 = 0. It easy to see that the 
duality principle is true indipendently by our ability to do the computations 
of the equations one gets from the perturbation series. 

In turn, the existence of a duality principle in perturbation theory means 
that is possible a perturbative analysis in different regions of the parameter 
space of the given equation. This situation could turn out to be very useful 
in quantum mechanics if one is able to obtain a dual Dyson series. This is 
indeed the case. 

So, let us consider the time-dependent Schrodinger equation 




being H (t) the Hamiltonian and h = 1 here and in the following. The Dyson 
series is a perturbative solution of this equation given by 

\*l>(t)) = (l ~ i I dtiH{tt) - f dh I" dt 2 H{t x )H(t 2 ) + ■■■) \ip(t )) 

V Jto Jto Jt J 

or, by introducing the time ordering operator T, 

|V(i))=Texp {-i jTV#(0) m )). 

The dual series can be obtained, through the duality principle, by assum- 
ing that the Hamiltonian H{t) has a discrete spectrum, that is H(t)\n,t) = 
E n (t)\n,t) with \n,t) the eigenstate corresponding to the eigenvalue E n (t). 
Then, the dual Dyson series is the one given in Ref. J10H , that is 

\t/j(t)) = U A (t)Texp (-i£dtH'(tf) |^(io)) 

being 

tt (+\ ijn.(t)-if dt'E n (t'), . 

U A {t) = 2_^e Jt o \n,t){n,t \ 

n 

the adiabatic unitary evolution operator, for the Berry phase j n (t) = 
{n, t\iMr\n, t) and 



H'(t) = - e- i ^( t )^W)e i ^^ ( ^ ( ^~^ W) (m > t|iR^|n J t)Kto)(n,*o|. 

This result proves that the well-known adiabatic approximation and its higher 
order corrections can be very effective in building asymptotic approximations 
to the solution of the Schrodinger equation, as is, on the other side, the Dyson 
series. 

Let us now consider a perturbed quantum system with Hamiltonian 

H = H + V(t) 

being H the Hamiltonian of the unperturbed system and V(t) the perturba- 
tion. In the interaction picture one has 

flj(t) = e iHot V{t)e- iHot . 

It is now possible to study the given system in different regions of the param- 
eter space through the Dyson series and its dual. In the former case we have 
standard textbook time-dependent perturbation theory. In the latter case we 
have to compute 

fli(t)|n,*>j = J BW(t)|n,t) J . 

But Hj(t) is just the interaction V(t) transformed by an unitary transfor- 
mation. Then, the eigenvalues En\t) are those of the perturbation V(t) and 
the eigenstates \n,t)i are just an unitary transformation away from the cor- 
responding eigenstates. These are the dressed states as generally computed 
in the current literature: It is just the leading order approximation of a dual 
Dyson series. But now we have a more general theory and higher order cor- 
rections can be computed. Beside, we realize why the dressed states are so 
effective in a strong field regime being obtained from the dual Dyson series 
that has a development parameter exactly inverse of the one of the Dyson 
series. 

It should be pointed out that both Dyson series and its dual can have 
the same kind of problems. One of the most important is surely the question 
of secularities: In any case, resummation of secular terms can be achieved 
through the renormalization group methods as pointed out, for quantum op- 
tics, in Ref. M. 

B. An example 

To give a clear insight of the working of the above analysis for a differential 
equation, let us consider the standard textbook example 

if>"(x) +a 2 (x)4>(x) = 0. 

that can be written in the form (the % factor is introduced just for convenience) 



dx ^ 0(x) J " ^ -ia 2 (x) J ^ cP(x) J ~ { > \ cf>(x) J ' { ' 

We can apply Dyson series and its dual. Dyson series is not normally 
applied to the above equation. Indeed, it gives the expansion 



6(x) 



J- a 



dX ' ( -ia\x') ) I X0 dX 'I X0 dX " a 2 (x') 



( <K*o) \ 



(13) 



In order to compute the dual Dyson series, we need to compute the eigen- 
vectors and eigenvalues of the matrix L(x). So, for the eigenvalue a(x) one 
has the eigenvector 

^ = y/-2ia(x) ( ~ia(x) ) (14) 
and for the eigenvalue —a(x) 

\2,x) = ^ (■),)■ (15) 
y/2ia(x) \ ia(x) J 

Then, one has for the Berry phases (2, x\i-^\2, x) = (l,x\i-4-\l,x) = and 
the unitary evolution operator (g) 

v , , _ 1 / «(aJo) cos (j^ dx'a(x')) sin (/* o dx'a(x')J \ 

y/a(x)a(xo) \ — a(x)a(x ) sin ( cZx'a(x')^ a(x) cos (j^ o dx'a(x')) i 

It is straightforward to see that 

gives the well-known Wentzel-Kramers-Brillouin-Jeffreys (WKBJ) result 

ijj(x) 1 cos f f dx'a(x')) -\ 2 sin f [ dx'a{x') \ . (18) 

y/a(x) \Jx ) \Ja{x) \Jx J 

In this derivation we have omitted the problem connected to turning points. 
We just note that, if there are points where a{x) = 0, Berry phases are no 
more zero as these are degeneracy points. 

This example shows the full power of the adiabatic approximation in find- 
ing asymptotic approximations to a given differential equation, without any 
requirement of slowly variation of the parameters of the equation. In the 
following we will show how to find higher order corrections too. 
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C. Duality and Berry's asymptotics 



Duality principle has been introduced in Ref. || to resolve problems both 
with infinitely small and large perturbations. As such, there is a region of 
the parameter space that is not possible to analyse by perturbation meth- 
ods. But, it is not difficult to realize that, as a by-product, an alternative 
solution to the Schrodinger equation for its unitary evolution through eq.(|5|) 
is obtained. This has no trivial consequences as, differently from the Dyson 
series, a superadiabatic scheme could be applied instead, as devised by Berry 



[11] that could give non-perturbative informations on the dual series. 

A superadiabatic scheme proves to be very useful when the full Hamil- 
tonian is considered with no a priori large or small parts, as shown in 
Ref. to describe stimulated Raman adiabatic passage by a three-level 
model. Indeed, the idea is to iterate the scheme to compute the adiabatic 
series giving Uj^if) and H'(t), by computing U' A {t) for H'(t), and the new 
Hamiltonian H"{t) through the eigenstates of H'{t). In principle, the pro- 
cedure can be repeated to the step one wants, giving the unitary evolution 
U{t) ~ U A {t)U' A {t)U%{t) ■ ■ ■ U^ n \t) and it is tempting to stop to a given step 
to obtain an approximation to the unitary evolution but, actually, the pro- 
cedure is shown to diverge. Anyhow, an optimal step n c exists for which 
an eigenstate basis set can be build by the approximated U{t) to approxi- 
mate the solution of the Schrodinger equation. Divergence is due to the fact 
that off-diagonal terms computed by the new Hamiltonians are systematically 
neglected. 

Indeed, to address the question of dressed states we consider a Hamiltonian 

like 

H = y <73 + V{t)a x 

being V(t) a generic perturbation, o\ and 03 Pauli matrices and ujq the level 
separation of the model. The regimes of interest are fully perturbative as 
V{t) is assumed to be very large. So, the initial Hamiltonian to apply the 
superadiabatic scheme is given, in interaction picture, by 

In this case, the superadiabatic scheme just stop to the second step as, at first 

step one has UA(t) = e % ~^ azt e~ trTl fa dt v ^ ) and, at the second step, U' A (t) = 
U A (t), so the product of unitary evolution operators is stopped and nothing 
new is obtained. Anyhow, the Berry's scheme can prove to be very useful 
in a non-perturbative regime, that is, when V(t) and ujq are of the same 
order of magnitude and exponentially small factors can be retained. Then, 
we can conclude that a superadiabatic scheme turns out to be useful in an 
intermediate regime, being in this way a bridge between the small and large 
perturbation theory linked in turn by the duality principle. This matter 
deserves further investigation. 



III. PERTURB ATIVE ANALYSIS OF THE JAYNES-CUMMINGS MODEL 



The Jaynes-Cummings model is widely used in quantum optics. Its Hamil- 
tonian, in the RWA, is given by M 



Hj C = uata + ^(|2)(2| - |1)(1|) + g(\2)(l\et + \l)(2\a) 

representing a two-level atom coupled with a single mode radiation of fre- 
quency <jj through the constant g. The reason to consider it here is that the 
exact solution is known and can be compared with the results of our pertur- 
bative analysis. 

In the interaction picture one has the Hamiltonian 



ffjg = 5 ( e ^|2)(l|a T +e- iAt |l)(2|a) 

being A = ujq — lo the detuning that here we assume different from for the 
sake of generality. As it can be seen from the form of H^jq, the critical param- 
eter in the model is the ratio ^. This means that an eventual perturbation 
series and its dual will have this parameter and its inverse as a development 
parameter. Now, we proceed to compute those series from the exact solution. 
The exact solution of the Schrodinger equation in interaction picture 



H 



Mi 



iJCWU-'-Qf 
can be found by looking for a solution in the form 

Wli = ci,n+i(t)\l,n + 1) + ca,„(t)|2, n) 



(21) 



(22) 



(23) 



(24) 



being n the photon number. So, the probability amplitudes are given by 

2igy/n + 1 



ci,n+l(t) = { Ci )Jl+ i(0) 



cos 



iA 



C2,n(i) = \ C2,n(0) 



COS 



iA 

a 



sin 



2ig\pn + 1 



-C2,n(0) sin 



Vt n t 



ci >n+ i(0) sin ( 



-iAt/2 



(25) 



being O n = yj A 2 + 7£ 2 and lZ n = 2g^J n + 1 the Rabi frequency. As expected, 
being A and g the only parameters, their ratio enters the only meaningful 
development parameter. The Dyson series is obtained by expanding the above 
solution in Taylor series of A = ^ giving till second order 



ci, n +l(t) = { ci, n+ i(0) 



l + »- 



. A 2 



(At + i(l 



X 



c 2 , n (0)(l-e- iA *) + O(A 



(26) 



C 2 ,n(i) = { C 2 ,n(0) 



1 - i^- (At + i(e iAt - 1 



- Cl , n+1 (0)(e lA *-l) + O(A 3 ) 



S 



(t) 



It is easy to see that at second order in the development parameter a secularity 
appears, that is a term that grows without bound in the limit t — > oo. In 
perturbation theory, unless we are not able to get rid of the secularity the series 
is not very useful. This can be accomplished through the renormalization 
group methods described in Ref. 0. But here, the problem can be easily 
traced back to the Taylor expansion of the functions sin(Vl + e 2 t) in e, having 
\/TTe2 = 1 + f + 0(e 4 ). So, we can eliminate it by simply substituting A 
with A + everywhere in the approximate solution into the exponentials of 
eq.©. 

It is not difficult to get back the result (|26| ) through the Dyson series (|3|). 
So, as expected, this series gives an analysis of the Jaynes-Cummings model 
when the detuning A is enough larger than the Rabi frequency TZ n . 

Now, let us repeat the above discussion in the opposite limit with the Rabi 
frequency larger than the detuning. Again, by Taylor expanding the exact 
solution one has 

Un.\ . i . /KnA 1 K n . (Tl ri 

1 —77 ism — 

2 J 2A 2 2 V 2 



c l,n+l(0) 
iC2,n(0) 



COS 



-t 



sm 



+ -r sin 
A 

2A^( Sin 



-t 



t cos 



+ 



(t) = C2,„(0) 



ic ln+ i(0) 



Tin 
COS j — t 

Sill j ——t 



i 



1 K n+ . (K n + 

— -77 rsm 1 

2A 2 2 V 2 



1 cos 1 

2 V 2 



+ o 



with the same problem of a secularity at second order. Indeed, this series 
can be obtained by the dual Dyson series (0) showing what could seem an 
unexpected result from the adiabatic approximation, but in agreement with 
the results of Ref. || . 

To compute the dual Dyson series we need the eigenstates and eigenval- 
ues of H^jq. It is easily found that for the eigenvalue gy/n + 1 we have the 
eigenstate 

\a,n,t) = -J=(e- iA *|l,n + l) + |2,n)) 
and for the eigenvalue —g\Jn + 1 we have the eigenstate 



\b,n,t) = -j=(\l,n+ 1) - e 



2,n» 



that are easily recognized as the dressed states of Ref. M for the Jaynes- 
Cummings model with a non-zero detuning. Berry phases are then easily 
computed to give 

7a = (a,n,t\i—\a,n,t) = — 



. , d . . A 
lb = {b,n,t\i—\b,n,t) = - — . 
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Then, after some algebra using the dressed states computed above, the unitary 
evolution operator (||) is given by, 



U (t) = e 



i £ t—igy/n+lt 



, n, t) (a, n, 0| + *+^v^+T* ^ ^ q ^ ^ Q | 



= cos (^r<J (e -< T*|l,n + l)(l,n + l| + e i T t \2,n)(2,n\) 
-isinf^^ (e~^*|l,n + l)(2,n\ +e i T*|2,n)(l,n + 1| 
that, for 1-0(0)) =c 1>n+1 (0)\l,n + l) + c 2>n (0)\2,n), gives 



IV-(*))j 



+ 



cos i ^t) c ln+ i(0) - ism (^t ) c 2 ,„(0) 



<<'s ( ^-t) c 2 ,„(0) - isin f-Tp* ) 'm./h-j 



e- l 2-*|l,n + l) 



that is the exact form of eqs.(|27|) when higher order terms beyond the leading 
one are neglected, i.e. when A — > oo, as expected from the results of Ref. |J. 

In order to go to higher orders, we have to compute H'(t) from eq.(|7|). 
Again, using the above expressions for the dressed states one gets 

H'{t) = -- [cos {TZ n t) (|1, n + n + 1| - |2, n)(2, n|) 

- isin(ft n i)(|l,ra + l)(2,n| - |2, n)(l, n + 1|)] 

so that, the first order correction to the leading order evolution operator Uo(t) 
of eq. (|32]) is given by 



Ux(t) = -%U (t) J dt x H' (h) = ijsin(^y?J (e^*|l, n + 1}(1, n + 1| -e^'l^raX^n] 

that gives the first order correction 

M(t))j =iisin^t) (e-^*ci, n+1 (0)|l,n + l) - e^c^^n)) 

again in agreement with the Taylor expansion as given in eqs.(p7[), to order 
t. So, in the same way we have at the second order 



U 2 (t) = -U (t) [ dtiH'ih) f h dt 2 H'(t 2 ) = 
J o J o 

i 2^ { [sm (^i) - ^ cos (^i) ] (e«#*|2, n> (1, n + 1| + e"^*] 1, n + 1) (2, n| 
- ^isin(-^i) (e- i ^|l,n + l)(l,n + l| +e4*|2,n)(2,n|)| 
then, one has 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 
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(e^* Cljn+1 (0)|2,n) + e- i $ t C2 tn (Q)\l,n + 1» 

(38) 

- i ^sin f^i) (e-^* Cl , n+1 (0)|l,n + 1) + &o^)% "»} ■ 

The agreement with the Taylor expansion as given in eqs.(^), to order is 
complete. 

As expected from results of Ref . |J , the adiabatic approximation and its 
higher order corrections turn out to be a nice method for asympotic analysis 
of the Schrodinger equation, being the dual of the well-known Dyson series 
and explaining in this way the nice working of the method of dressed states 
currently used in quantum optics. No slowly varying of the parameters of the 
Hamiltonian is involved as one could expect for the adiabatic approximation. 



IW))j 



'2A 2 



SID [ — t 



— tcosf— t 



IV. PERTURB ATIVE ANALYSIS OF MODELS FOR HIGH-ORDER HARMONIC 

GENERATION 

A. Models 

Several models are currently used to account for high-order harmonic gen- 
eration. The first model considered |2| has been a one-dimensional model 
described by the Schrodinger equation 



i_d 2 

2dx 



-^2 + V ( x ) ~ xe (t)smu; L t 



.dV(x,t) 
' dt 



being eo(t) a function taking in account the time to rise the laser field to 
its maximum value, lol the frequency of the laser field and V(x) a simple 
representative binding potential for the atom. A choice currently found in 
literature is V{x) = — ^ppp ■ Beside numerical methods that are very com- 
puter demanding, other methods as Floquet theory have also been applied 
[13] for the full three-dimensional case. A fruitful understanding of harmonic 



generation through semiclassical ideas has also been yielded in Ref. [14]. By 
these semiclassical results, a non-perturbative quantum model has been ob- 
tained [[lq ]. Beside, a first approach by second quantization has also been 
given where a hint was put forward that harmonic generation is a first order 
effect [16[. Analytical expression are barely given as all these models have 
been solved or numerically either non-perturbatively so to require at some 
step numerical computation. Another model is a simpler two-level system 
described by the hamiltonian 



H 



u 



(|2)(2|-|l)(l|)-xe (t) 



sin u^t 

COS LOjjt 



(39) 



(40) 



and 



x = -d 12 (|l)(2| + |2)(l|) 



(41) 
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being d\2 the matrix element of the atomic dipole. This model is well-known 
in quantum mechanics. A first hint to a strong coupling perturbative solution 
was given by Meystre § that used it to describe an atom inside a Fabry- Perot 
cavity. The series till first order and the way to compute higher orders for 
strong coupling were finally obtained in Ref. |7j where it was shown that a 
shift of the levels of the atom occurs. 

Indeed, this two-level model seems very effective in describing high-order 
harmonic generation too. The two physical situations of a Fabry-Perot cavity 
strongly coupled with an atom and an atom in a strong laser field seems de- 
scribed by the same hamiltonian. But this should not come out as a surprise. 
What really matters here is the existence of the shift of the energy levels of 
the atom in these situations that, for the case of high-order harmonic genera- 
tion can change the spectrum of hyper-Raman lines and so, can be measured 
experimentally. 

Beside, as we are going to show, the leading order solution found by 
Meystre is not enough to get the power spectrum computed through the 
Fourier transform of the equation 

x(t) = (V(t)\x\V(t)) 

In fact, by the dual Dyson series one can see that high-order harmonic gen- 
eration is actually a first order effect. In this way we are able to reproduce 
the results obtained in Ref. [|| by Floquet method, but having an analyti- 
cal expression to be compared with experiments. As a by-product we have 
that the hyper-Raman lines can be shifted. Through this approach the com- 
putation can be pushed to any order, coping always with definite analytical 
expressions. 

The model (|39| ) can also be treated by this approach. Indeed, an applica- 
tion to multiphoton ionization has been found by Salamin [j^j. The leading 
order solution should be written as 



tl>(x,t)ne i *fi deeo W'* nUl '*<l> n (x) 



being 



id 2 



4> n (x) = E n (j) n (x) 



It easy to see that probability transitions given by w mn (t) = 
/-to dx())m(x)ip(x,t) are not trivial and can be computed also for the con- 
tinuos part of the spectrum. But, as we are going to show using the two-level 
model and as can be seen by the look of the leading order solution (|43|), we 
need to compute the first order correction to it to account for high-order har- 
monic generation. We do not pursue the study of this model further here, 
as the two-level model can give a satisfactory account of all this matter in a 
simpler way. We just note that in this way, more complex models than that 
of eq. (|39|) , through perturbation methods, could be taken into account. 
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B. Perturbative analysis for high-order harmonic generation 



To fix the ideas, we consider the two-level model of Ref. ||, that is, eq.([40|) 
with a cosine perturbation. Dyson series using probability amplitudes and its 
dual solution to first order of this model through operatorial methods were 
given in Ref. |7j. So, we avoid the analysis by the Dyson series of this model 
discussed in depth in and Refs. therein. Instead, we use the dual Dyson 
series to show that it is equivalent to the operatorial method used in |7| and 
presented initially in Ref. |l8| . 

The rising of the laser field accounted for by the function eo(t) is taken 
as istantaneous to make the computations simpler, that is, we take eo(i) = 
=constant. 

In interaction picture, the Hamiltonian ( f40|) is given by 



Hi = fld 12 cos uj L t{e~ luJot \l)(2\ +e^°'|2)(l|) (45) 

Then, computing the dual Dyson series, for the eigenvalue Qdi2 cos uit we 
get the eigenvector 

|M> = ^(e W |2) + |l)) (46) 

and for the eigenvalue — Q-d\2 cos oj Lt the eigenvector 

|o J t> = -L(|2)-e-*-»*|l». (47) 

These are the dressed states for this model. The corresponding Berry phases 
are given by 

76(t) = y 



Ja(t) = - — ■ 



(48) 



It is interesting to note here, that Berry phases originate from the energies of 
the levels of the unpertubed atom. 
All this gives the unitary evolution 

U Q (t) = e-^ t e t9 ^T sinuJLt \a,t)(a,0\ + e ^V^ sin ^|6, 0| (49) 
that yields in terms of the bare states |1) and |2) 

Uo (t) = cos ^sin^i) (e-^|l)(l|+e^|2)(2|)-zsin (^riawrf) ( e -^|l)(2| + e^|2)(l|). 



(50) 



We can reformulate the above operator as a matrix by taking for the bare 
states 
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so to have 
U (t) = 



e *^t cos f£Mui sin^ -ie*^* sin sin w L t 
-ie - *^* sin sina; L t) e - *^* cos sin w L t 



It is not difficult to see that the above operator can be rewritten through the 
Pauli matrices <7i,<72>0'3 as 

then, by eliminating the prefactor due to interaction picture, we are left with 
the leading order result of Ref. |?| for the wave function 

\*(t))^e~ tCT1 ^r s ' muJLt \*(o)). 

In the same way, we can compute higher order corrections to the above by 
computing H'(t) for the dual Dyson series. In the bare states, using again 
the dressed ones, one has 



2 



cos j 2^^-smu L t\ (|2)(2| - |1)(1|) -isin (^^-smuj L t ) (|2)(l. - 1)(2|) 



that is, 



if (t) = ye <"£ a 3 e w £ 



in agreement with the computation of the first order correction computed 
through operatorial methods in Ref. 0. The two series are identical as it 
should be expected. 

So, the solution for the system ( |40| ) till first order can be written as 

It/.W — iai — —smu)rt 

|*(t)) = e <"£ x 

[ I ~ i-^-Jo ta 3 



2 J \ LU L 
^ /2f2di2\ sin(2ncj£i) 

iu 2^ ^2n ( — ] 0"3 



n=l 
oo 



\ uol J 2nu>L 



. ^ T /2ndi 9 \ cos((2n + lWt) - 1 

+ ZW 2^ J 2 n+1 pr —TT (T2 

^o V / (2n + l)w L 

+ HI*(o)>. 

where use has been made of the operatorial identity 

oo oo 

e ±i<Tfe2sin * = J (z) + 2 ^ J 2n (z) cos(2n</>) ± 2ia k ]T J 2n+ i(z) sin((2n + 1)0) 

n=l n=0 
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with <7fc one of the Pauli matrices and J n (z) Bessel functions of integer order. 
The secular term in eq.(|57|) can be resummed away by renormalization group 
methods, as shown in Ref. JiJ, giving the renormalized levels of the atom 
in the laser field. Then, the solution one has to use to compute the power 
spectrum is 

I t /,\\ —MTl — SinoJrt 

| = e u l x 

I - iujq > J 2n — 03 

^ \ cjl J 2nuj L 

. ^ T f2Qd 12 \ cos((2n + 1W) - 1 

+ IUJQ J 2n+1 77, — TT 0"2 

j^o \ u L J (2n + l)w L 



+ ...] e 2 "V ^ 7 |^(0)). 

It easy to see that if we just limit our analysis to eq.(^), we are not able 
to obtain the spectrum of the harmonics. In fact, one would have from eq. (fl2"|) 
x(t) = (*(t)|a;|*(t)) = — di2 (°) I °"l I * (0) ) =constant. Instead, using eq. (||) 
one has at first order 



x(t) = -d 
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c 2 c\e- iuJ0Rt + c* 2 c ie iuJ0Rt 



2 i \2\ t f2Qd 12 \ cos((2n + l)u L t) - 1 

+ ci - c 2 )u 2^ J 2n+i : — ' Y\ 

^ V w L J (n + ±W 



n=0 v ~ u ' y t 2 

oo 



+ i(c* 2Cl e iuJ0Rt - c 2 c\e- iulORt )uJo ^ J- 



In 

n=l 



/2r2d 12 \ sin(2no;it) 



V lol J nu) L 



being loqr = luoJo ^ 2 ^ 12 j the renormalized separation of the two levels of 
the atom and introducing the population distribution c\ and c 2 for the bare 
levels of the atom through the initial state |^(0)). This is exactly the form 
one must have for the high-order harmonic generation using this model, as 
showed in Ref. || by the Floquet method. In fact, we have odd harmonics of 
intensity (|ci| 2 — |c 2 | 2 ) 2 , while the latter term is due to the hyper-Raman lines 
at uoii±2nuiL of intensity |c 2 | 2 |ci| 2 . But now, an explicit analytic expression 
for the power spectrum is given, so that, a clear understanding of all the 
parameters involved into atom and strong laser field interaction is obtained. 
Particularly, we have an exact expression for the shift of the levels of the atom 
that now we know how to measure: If we take initially |ci| = |c 2 | we will be 
just left with hyper-Raman lines into the spectrum. These lines can be shifted 
by varying the ratio 2 ^ 12 , given by the parameters to be controlled into the 
experiment. By observing how these lines move we can obtain a measure of 
the shifts of the atom levels. By the expression of x(t), it is clear that the 
only lines that can be moved are indeed the hyper-Raman lines. 

The advantges above the Floquet method used in Ref. || are evident as 
we have an explicit analytical formula for the spectrum with all the functional 
dependencies on the parameters entering into the model explicitly expressed. 
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V. CONCLUSIONS 



The way dressed states are currently treated in quantum optics becomes 
to compute the eigenstates and eigenvalues pertaining or the interaction term 
either the full Hamiltonian of a given system. In this paper we have shown how 
the dual Dyson series can give both an understanding of this approach and a 
tool to improve it. Indeed, the existence of a dual Dyson series permits, as we 
have shown, a perturbative study of atom-field interaction in different regions 
of the parameter space of a given model. The case of the Jaynes-Cummings 
model, fully exploited in this paper, is examplary in this sense. Beside, the 
dual Dyson series is nothing else than the quantum adiabatic approximation 
and its higher order corrections that in this way prove to be a very powerful 
tool to obtain asymptotic approximations to the Schrodinger equation. 

The theory is applied to the analysis of models for high-order harmonic 
generation giving an explicit expression for the power spectrum. In this way 
new experiments can be thought where, properly changing the parameters 
involved, measures of the energy levels of the atom in interaction in a strong 
laser field could be accomplished. Beside, a control on the amplitudes of the 
harmonics and the position of the cut-off on the spectrum could be obtained. 
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